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ABSTRACT 

The optical afterglow spectrum of GRB 050401 (at z = 2.8992±0.0004) shows 
the presence of a damped Lya absorber (DLA), with log N H i = 22.6 ± 0.3. This 
is the highest column density ever observed in a DLA, and is about five times 
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larger than the strongest DLA detected so far in any QSO spectrum. From the 
optical spectrum, we also find a very large Zn column density, allowing us to 
infer an abundance of [Zn/H] = —1.0 ± 0.4. These large columns are supported 
by the early X-ray spectrum from Swift-XRT which shows a column density 
m excess of Galactic) of logiV H = 22.21+°;°^ assuming solar abundances (at 
z = 2.9). The comparison of this X-ray column density, which is dominated 
by absorption due to a-chain elements, and the H I column density derived 
from the Lja absorption line, allows us to derive a metallicity for the absorb- 
ing matter of [a/H] = —0.4 ± 0.3. The optical spectrum is also substantially 
reddened and can be well reproduced with a power-law with SMC extinction, 
where Ay = 0.62 ±0.06. But the total optical extinction can also be constrained 
in a way which is independent of the shape of the extinction curve: from the 
optical-to- X-ray spectral energy distribution we find, 0.5 < Ay < 4.5. However, 
even this upper limit, independent of the shape of the extinction curve, is still 
well below the dust column that is inferred from the X-ray column density, i.e. 
Ay = 9.1^^5. This discrepancy might be explained by a small dust content with 
high metallicity (low dust-to-metals ratio). 'Grey' extinction cannot explain the 
discrepancy since we are comparing the metallicity to a measurement of the total 
extinction (without reference to the reddening). Little dust with high metallicity 
may be produced by sublimation of dust grains or may naturally exist in systems 
younger than a few hundred Myr. From these results it is clear that dust ex- 
tinction properties in GRBs derived from comparisons of optical reddening and 
metallicity are unreliable. 

Subject headings: gamma rays: bursts — X-rays: general — galaxies: ISM — galaxies: 
high redshift — quasars: absorption lines — ISM: dust, extinction 

1. Introduction 

The largest H I column density absorption line systems [iV(H I) > 2 x 10 20 cm -2 ], known 
as damped Lya absorbers (DLAs), were first observed in QSO spectra, and originate in the 
H I regions of intervening galaxies (see Wolfe, Gawiser, & Prochaska 2005, for a recent 
review). Because background QSOs probe random sightlines through intervening galaxies, 
QSO-DLA systems can in principle be used as an unbiased tracer of the neutral gas at high 
redshift, gas that provides the fuel for star formation at these epochs (Lanzetta, Wolfe, & 
Turnshek 1995; Wolfe et al. 1995), though complications regarding a bias against dusty lines 
of sight may possibly exist (Vladilo & Peroux 2005; Akerman et al. 2005). 
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Typically, metallicities are low in QSO-DLAs compared to solar metallicity, ([X/H] ~ 
—1.5 at z ~ 3) with a large scatter (Kulkarni & Fall 2002; Prochaska et al. 2003). 

DLAs are also common in long duration 7-ray burst (GRB) afterglow spectra (Jensen 
et al. 2001; Hjorth et al. 2003a; Jakobsson et al. 2004b; Vreeswijk et al. 2004), typically with 
much higher column densities than observed in QSO-DLAs (Vreeswijk et al. 2004). It seems 
highly plausible, given that long-duration GRBs are known to have massive star progenitors 
(Stanek et al. 2003; Hjorth et al. 2003b), that the extremely large neutral hydrogen column 
densities observed in many GRB-DLAs are related to the GRBs' star-forming regions (e.g. 
Bloom, Kulkarni, & Djorgovski 2002). This means that while the highest redshift absorption 
system in an afterglow spectrum belongs to the GRB host, and is therefore unlikely to be 
useful as a random probe of intervening matter, it can instead provide rich information 
on the sites of active star-formation in the high-redshift universe (the mean and median 
redshifts of recent samples of GRBs are now very high, around z = 2.8, and the first GRB 
at z > 6 has recently been found; Jakobsson et al. 2006; Kawai et al. 2006; Tagliaferri et al. 
2005; Price et al. 2005; Haislip et al. 2006; Watson et al. 2006). At the same time, afterglow 
spectra can still be used as random probes for intervening absorbers at lower redshifts in 
the same way as QSOs. The largest Lja absorption columns in DLAs are found in GRB 
afterglows: GRB 030323 with logN H i = 21.90±0.07 (Vreeswijk et al. 2004) and more recently 
logNm = 22.1 ±0.1 in the afterglow of GRB 050730 (Starling et al. 2005; Chen et al. 2005). 

Further data on the birth matrix of GRBs are provided by studies of the afterglow's 
optical extinction and soft X-ray absorption. The optical studies provide an opportunity to 
constrain the extinction law for star- forming galaxies at high redshift, though the work may 
be complicated by the possibility of destruction of dust grains by the GRB. Indeed it has 
been proposed that GRB afterglows are subject to a 'flat' or 'grey' extinction law due to a 
change in the grain size distribution caused by the sublimation and breaking of dust grains 
(Waxman & Draine 2000; Fruchter, Krolik, & Rhoads 2001; Perna, Lazzati, & Fiore 2003), 
i.e. where there is evidence for little spectral reddening but where the presence of a significant 
dust column is inferred. Proposed evidence for this includes large metal column densities 
observed in afterglows which exhibit little curvature in their optical afterglow spectra. The 
metal column densities are derived either from optical absorption lines (where dust depletion 
effects have also been observed, Savaglio, Fall, & Fiore 2003; Savaglio & Fall 2004), or from 
the soft X-ray absorption (Stratta et al. 2005, 2004), where the most important contribution 
comes from oxygen K-shell absorption. This apparent discrepancy between the metal column 
density and the optical reddening was first observed in a careful and seminal paper by Galama 
& Wijers (2001), using very low quality X-ray spectra from BeppoSAX. 

The only alternative to a flat extinction law appears to be a non-universal dust-to- 
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metals ratio. This would explain a high metal column density with little spectral reddening 
without reference to non-standard extinction laws, since there is simply less dust. A low 
dust-to-metals ratio appears a very natural possibility for GRB host galaxies, given their 
high redshifts (Jakobsson et al. 2006) and very young stellar populations (typically less than 
two hundred Myr, Le Floc'h et al. 2003; Christensen et al. 2004; Gorosabel et al. 2003, 2005). 
This is because metal enrichment of the interstellar medium depends on supernovae (SNe) 
which have short lifetimes, while the overwhelming mass of dust (at least in the current 
epoch) is not produced by SNe, but by asymptotic giant branch (AGB) stars (Andersen, 
Hofner, & Gautschy-Loidl 2003) that have lifetimes of at least one hundred Myr before even 
the most massive stars begin to produce quantities of dust (Ferrarotti & Gail 2005). 

To distinguish between these possibilities, we require a direct measure of the hydrogen 
Lya, soft X-ray, and optical metal absorptions in a single GRB afterglow, observations that 
have been elusive because of the conflicting requirements of a relatively low redshift for an 
accurate determination of the intrinsic soft X-ray absorption, and high redshift to move 
hydrogen Lja to observable wavelengths. 

In this paper we present X-ray and optical observations of GRB 050401 (§ 2), the result- 
ing spectra and lightcurves (§3), and the implications of these observations (§4). 

A cosmology with f2 m = 0.3, = 0.7 and H = 75kms _1 Mpc -1 is assumed through- 
out. Error ranges quoted are 68% confidence intervals for one parameter of interest, unless 
otherwise indicated. 
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Table 1: Optical and near-IR imaging observations 

Date (UT) Telescope Mag Time since trigger (ks) 



R-band 



April 


1.9880 


MAO 


22. 


.7(2) 


33.75 


April 


2.0720 


TNG 


23 


.0(1) 


41.01 


April 


2.2100 


VLT 


23 


,27(9) 


52.93 


April 


2.2243 


NOT 


23 


.0(2) 


54.16 


April 


2.2563 


VLT 


23 


,31(8) 


56.93 


April 


2.3004 


VLT 


23 


.44(8) 


60.74 


April 


2.3183 


D1.5 


23 


.5(2) 


62.29 


April 


2.3438 


VLT 


23 


.5(9) 


64.49 


April 


3.2414 


NOT 


23 


,9(2) 


142.0 


April 


8.2116 


NOT 


25. 


.3(4) 


571.5 


April 


14.5947 


Subaru 


26. 


.0(1) 


1123 



J-band 

April 1.6402 UKIRT 18.6(2) 3.698 
H-band 

April 1.6346 UKIRT 17.9(1) 3.214 

April 1.6457 UKIRT 18.1(2) 4.173 

April 1.6513 UKIRT 18.0(1) 4.657 

April 2.3521 CTIO 20.3(1) 65.21 

K-band 

April 1.6242 UKIRT 16.88(8) 2.316 

April 1.6284 UKIRT 17.2(1) 2.678 

April 14.2742 VLT > 21.5 1095 



Note. — Uncertainties in the last significant digit are listed in parentheses after that digit. Observations of 
the afterglow of GRB 050401 were obtained at the ESO Very Large Telescope (VLT), the 8.2 m Subaru tele- 
scope, the 2.5m Nordic Optical Telescope (NOT), the Danish 1.54m telescope (D1.5), Telescopio Nazionale 
Galileo (TNG), the 1.5m telescope at the Maidanak Astronomical Observatory in Uzbekistan (MAO), the 
Blanco 4.0m telescope at the Cerro Tololo Inter- American Observatory (CTIO), and the 3.8m UK Infrared 
Telescope (UKIRT). 
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2. Observations and data reduction 

On 1 April 2005 at 14:20:15 UT Swift' 's Burst Alert Telescope (BAT) triggered on a 
multi-peaked burst with a duration of 33 ± 2s (Barbier et al. 2005; Sakamoto et al. 2005). 
A rapid autonomous slew of the satellite resulted in the X-Ray Telescope (XRT) being on 
source within about two minutes and a bright, fading source was discovered in the BAT 
error-circle (Angelini et al. 2005). XRT revisited the source again on 6, 7 and 8 April. 
We acquired the XRT data from the archive and used the level 2 events files to construct 
the lightcurves and spectra, using a 40 pixel wide source extraction box for the Windowed 
Timing (WT) mode data with a background region box 4' away. A 12 pixel radius circle was 
used to extract the source counts in Photon Counting (PC) mode; an annulus around the 
source region was used to extract a background in this case. Data in PC mode with a source 
count rate above 0.4 counts s -1 are likely to be affected by pile-up and were not used, as good 
statistics on the bright afterglow were available from the WT mode data. HEASOFT 6.0 was 
used to analyse and reduce the data with the latest calibration files. Ancillary response files 
were produced with xrtmkarf using 'inarflile=CALDB', but the fit to the data at energies 
close to 0.5 keV was poor. Running xrtmkarf with 'inarflile=NONE' produces an ancillary 
response file without using a modified on- axis ancillary reponse as input, instead using the 
mirror on-axis effective area and the filter transmission file directly. The ancillary response 
produced in this way improved the fit at low energies, though the results presented here 
are not strongly dependent on the choice of response file, in particular a very large excess 
column density is detected regardless of which ancillary response is used. 

Optical and near-IR observations of the optical afterglow of GRB 050401 (McNaught & 
Price 2005) were secured at several telescopes (see Table 1). Starting at 5:16 UT on April 
2, 2005 (14.9 hours after the burst trigger), 8 spectra of 1450s each were acquired with the 
FOcal Reducer/low dispersion Spectrograph 2 (FORS2) instrument on Antu of the Very 
Large Telescope (VLT) at Cerro Paranal in Chile, using the grism 300V and order sorting 
filter GG375, and a slit width of 1". This setting results in a binned pixel scale of 3.24 A, and 
a FWHM resolution of 9.5 A at 4000 A and 11 A at 9000 A. After cosmic ray removal using 
the van Dokkum (2001) algorithm, the spectra were reduced in the standard fashion, using 
IRAF. As the night was photometric, the GRB afterglow spectra were flux calibrated using 
the spectrophotometric standard star G138— 31. Finally, the spectra were corrected for the 
Galactic foreground extinction estimate of E(B— V)= 0.065 mag of Schlegel, Finkbeiner, & 
Davis (1998). The resulting normalized spectrum and its la error spectrum are shown in 
Fig. 1. 
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Fig. 1. — Lower panel: VLT/FORS2 spectrum of GRB 050401 obtained 15 hours after the 
burst when the afterglow magnitude was R = 23.3. There are two absorption systems at 
z ~ 2.9 (black labels, down arrows) and z ~ 2.5 (grey labels, up arrows). The DLA at 
z ~ 2.9 is the strongest ever observed and has a metallicity at least ~ 0.1 times the solar 
value ([Zn/H]= —1.0 ±0.4). The la error spectrum is also shown. The model superimposed 
on the observed spectrum corresponds to a fitted column density of logiV(Hl) = 22.6 ± 0.3 
together with possible H 2 Lyman-band lines from the J = and 1 rotational levels, with 
logA^(H 2 ) = 21 in both J. Upper panel: The 2 dimensional spectrum covering the same 
wavelength range as the 1 dimensional spectrum in the lower panel. The spectrum in the 
middle is the afterglow. Above and below are the spectra of 2 other objects on the slit. The 
spectrum of the object above the afterglow spectrum extends far bluewards of where the 
afterglow spectrum falls off. The spectrum below the afterglow spectrum is most likely due 
to a z = 2.65 Lya emitter. 
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3. Results 

3.1. The X-ray afterglow 

The WT-mode spectrum (Fig. 2) was fit with a power-law with fixed Galactic absorption 
(4.8 x 10 20 cm~ 2 , using the nh ftool, Dickey & Lockman 1990). The fit using this model was 
poor (x 2 /dof = 629/263), leaving large negative residuals at low energies, and so an extra 
parameter, additional absorption at z — 2.8992 (the redshift of the burst — see below), was 
added to the fit. This improved the fit significantly (x 2 /dof = 297/262), and gave best-fit 
values of photon spectral index T = 1.89 ± 0.03 and equivalent hydrogen column density 
TVh = 1.6 ± 0.1 x 10 22 cm~ 2 (log Nn = 22.20 ±0.03). Using the ancillary response file created 
with 'xrtmkarf inarfnle=CALDB', the best-fit values were similar — r = 1.85 ± 0.03 and 
N H = 1.6 ± 0.2 x 10 22 cm~ 2 . The absorbing column density is fit using the absorption cross- 
sections of Morrison & McCammon (1983), which assumes solar abundances. Splitting the 
WT-mode data in three equal duration bins revealed no significant change in the spectrum 
(absorbing column or spectral slope). Also the PC- mode data, obtained several hours after 
the burst, are consistent with the same absorbing column and slope (Fig. 2). An analysis 
of the first minute and the first five minutes of WT-mode data revealed no variation in the 
absorbing column density greater than the la level. The absolute systematic uncertainty in 
the flux normalisation of the spectrum is small, of the order of a few percent. 



-10- 




Energy (keV) 

Fig. 2. — X-ray spectrum of the afterglow of GRB 050401 folded through the instrument 
response (WT mode — crosses, PC mode — triangles). The best-fit model is plotted as a 
solid line. The model used is a power-law with absorption fixed at the Galactic level, plus 
absorption at z — 2.8992 (Nu = 1.6 x 10 22 cm -2 at solar abundances). Residuals from this 
model are plotted in the lower panel in units of standard deviations of the data. 
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The level of the Galactic foreground absorption may be uncertain up to ~ 20% (Dickey 
& Lockman 1990). Including this uncertainty leads to a more conservative error estimate of 
A^h = 1-6 ± 0.3 x 10 22 cm -2 for the equivalent hydrogen column density. 

Although the X-ray column density is usually expressed in terms of the equivalent 
hydrogen column density, this assumes solar abundances; the soft X-ray absorption is in 
fact a direct measure of the total metal column density in the solid and gas phases. The 
biggest contribution to the absorption comes from the oxygen K-shell absorption edges, with 
O I at 0.52 keV. There is also a significant contribution from other a-chain elements. For 
instance, at z ~ 3 the oxygen and other a-elements are responsible for two-thirds of the 
absorption in a typical X-ray CCD spectrum (e.g. Swift-XKT) assuming solar abundances. 
The X-ray column density is therefore effectively a measure of the a-element column (see 
also Morrison & McCammon 1983; Galama & Wijers 2001; Turnshek et al. 2003). It is very 
often impossible to determine the ionisation state of the X-ray absorber due to 1) Galactic 
absorption, 2) the fact that at moderate redshifts the oxygen edge is usually shifted out 
of the bandpass of most X-ray detectors, and 3) the low resolution of CCD X-ray spectra. 
This, combined with the fact that the difference in the measured absorbing column density 
is expected to be small between neutral and moderate ionisation states means that column 
density measures from soft X-rays are given for a neutral column. The fact that we do not 
know the ionisation state of the X-ray absorber leads, at worst, to a fairly small underestimate 
of the metal column. In this case there is evidence of a very large, neutral column density 
from the optical spectrum (see below), which makes it plausible to assume that the fraction 
of fully ionized gas is negligible. 

The 0.2-10.0 keV lightcurve (Fig. 3) is not well fit by a single power-law (x 2 /dof = 
983.5/65) but a broken power-law improves the fit significantly (x 2 /dof = 81.6/63) with a 
break at 4090 ± 350 s, from a decay index a X i = 0.58 ± 0.01 to a X2 = l-39±g;g|. There is 
no significant change in the spectrum from one side of the break to the other, although the 
spectrum on the late side of the break has low number counts and is not as well-constrained 
as the early spectrum. 
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Days since burst 

Fig. 3. — X-ray and optical lightcurves of the afterglow of GRB 050401. The R-band 
lightcurve (left axis, diamonds) is well fit by a single power-law with index ao — 0.82. 
This fit is consistent with data beginning at less than one minute after the burst, from 
Rykoff et al. (2005, open triangle) and McNaught & Price (2005, filled triangle) though they 
are not included in the fit. The observed X-ray lightcurve (right axis, WT mode — crosses, 
PC mode — boxes) requires a break at 4090 s with a x changing from 0.58 to 1.39. The epoch 
to which the data were interpolated to produce the SED is indicated with a dashed line. 
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3.2. Optical Photometry and Spectroscopy 

The optical lightcurve (Fig. 3) is well fit with a single power-law (x 2 /dof = 11.5/11) with 
optical decay index ao — 0.82 ±0.03. This fit is consistent with data from less than a minute 
(the burst was observed very rapidly by ROTSE-IIIa, Rykoff et al. 2005) to nearly two weeks 
post-burst. Such a uniform decay is unusual for optical afterglows that typically display one 
or more breaks in this time interval. Furthermore, this optical decay slope is significantly 
different from both X-ray decay slopes. The uniformity of the decay, however, and the 
similarity with the H-band decay (ojh = 0.76 ± 0.04) means that we can have confidence in 
the correction of the optical/NIR data to an epoch common to the X-ray observations. This 
confidence is important in order to reconstruct the spectral energy distribution (SED) of the 
afterglow correctly. 

A striking feature in the Galactic-extinction-corrected, flux-calibrated optical spectrum 
is that it is very red compared to typical optical afterglow spectra (Fig. 4). The best-fit 
spectral index {F v oc u" 13 ) is (5q = 2.8 (fitting the spectral range 5800-9000 A). This fact, 
coupled with an optical-to-X-ray spectral index (/?ox) of 0.2, strongly suggests that the 
optical spectrum is reddened by dust. SMC, LMC and Galactic extinction laws (Pei 1992) 
were fit to the optical spectrum assuming an underlying power-law continuum with the 
normalisation, the underlying power-law and the extinction as free parameters. As with all 
GRB afterglows studied to date (but see Vreeswijk et al. 2005), the lack of a 2174 A bump 
in the spectrum means that of these three, the SMC extinction law provides the best fit 
(e.g. Vreeswijk et al. 2004; Jakobsson et al. 2003; Savaglio & Fall 2004), yielding an intrinsic 
power-law spectral index, j3o = 0.5±0.2 and Ay = 0.62±0.06. We do not of course constrain 
the shape of the extinction curve very well, so that this measure is heavily dependent on the 
SMC curve. 
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Fig. 4. — The spectral energy distribution of the afterglow of GRB 050401. The best-fit 
single power-law (f3 = 0.89) to the X-ray data is plotted and extrapolated to optical-NIR 
wavelengths. All data has been interpolated to a common epoch of 17.9 minutes after the 
burst using the R-band lightcurve. There is no evidence for a break in the X-ray power-law 
spectrum. To derive the minimimum extinction required to dim the extrapolated X-ray 
power-law to the fluxes observed in the optical-NIR, the most conservative likely broken 
power-law (/3 = 0.39) is also plotted, with the break set at as high an energy as possible, 
at 0.4 keV (setting the break at lower energies requires larger extinction, while a higher 
break energy is ruled out by the X-ray spectrum). Fitting SMC extinction to the optical 
continuum spectrum (rebinned here by a factor of 10) assuming an underlying /3 = 0.39 
broken power-law, gives a good fit with Ay = 0.67 and is plotted as a dashed line. The SMC 
extinction that would be required to account for the R-band photometric data assuming the 
single power-law (f3 = 0.89) continuum extrapolated from the X-rays, is A v = 1.45, and 
this is plotted as a dotted line for comparison. However this is clearly not compatible with 
the observed optical spectral shape. Inset. The maximum and minimum extinctions as a 
function of wavelength. These limits for the total extinction are derived from the single and 
broken power-laws (f3 = 0.89 and f3 = 0.39 outlined above) divided by the observed data. 
Circles represent the observed NIR data, triangles the observed R-band, and the solid curve 
the continuum from the spectral observations. By contrast, the extinction expected from 
the soft X-ray absorption measurement is plotted as a filled square. The inconsistency of 
the optical/NIR limits with the extinction expected from the X-ray absorption shows that 
there is no room in the SED to allow for a total extinction that matches the high metallicity, 
suggesting a low dust-to-metals ratio. 
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We used this best-fit reddened power law and its extrapolation at wavelengths shorter 
than 5800 A to normalise the optical spectrum. A striking feature in the spectrum is the 
very large absorption edge blueward of about 5400 A (see Fig. 1, upper panel). Given the 
detection of two metal absorption line systems in the spectrum with redshifts of z ~ 2.5 
and z ~ 2.9 (see Table 2) we interpret the observed very large absorption trough as the 
signature of a very strong H I Lya absorber at z ~ 2.9. The correctness of this identification 
is based on the observation of the red damped wing of the Lya line all the way from 4800 to 
5300 A (see Fig. 1, lower panel). Fitting a damped Lya profile to this damped wing yields 
a neutral hydrogen column density of logiV(Hl) = 22.6 ± 0.3. We adopt a conservative la 
error of 0.3 dex which takes into account the main source of uncertainty which is related 
to the continuum placement, not to the RMS error from fitting the Voigt profile. The 
measurement of the rest-frame equivalent width on the red-half side of the line (removing 
the contributions of narrow absorption lines), ~ 75 A, independently confirms the fit result. 
This makes this system the largest DLA ever observed, larger even than that observed in 
GRB 050730 (Starling et al. 2005; Chen et al. 2005) and GRB 030323 (Vreeswijk et al. 2004), 
and a factor of five larger than observed in any QSO-DLA 1 . Given the very large H I column 
density, we checked for the presence of H 2 absorption lines at the GRB redshift. In Fig. 1 
(lower panel), we display a model of the Lyman-band lines of H 2 from the J = and 1 
rotational levels, assuming logA^(H 2 ) = 21 in both J. This would correspond to a molecular 
fraction of log / ~ -1 with / = 2iV(H 2 )/(2iV(H2) + iV(Hl)). However, we cannot claim 
a detection given the poor quality of the afterglow spectrum bluewards of 4500 A. Among 
the observed metal absorption lines, the Zn II doublet at z^ = 2.8992 is unambiguously 
identified. The observed equivalent width of the Zn II A2026 line is 5.5 ± 0.9 A. Since this 
line may be saturated, the corresponding Zn metallicity, [Zn/H]= —1.4, derived from the 
optically thin limit approximation, should be considered as a lower limit. We tentatively 
performed simultaneous Voigt-profile fitting of the Zn II doublet lines together with carefully 
selected, weaker and unblended lines from Cr II, Si II and Fe II (see Fig. 5). Assuming a 
single component and the same (turbulent) broadening parameter value for all ions, we find 
a best-fit redshift of 2.8992 ± 0.0003 and broadening parameter b = 67 ± 35kms _1 . The 
best-fit column densities are given in Table 3. Including the Mg I A2026 line in the fitting 
process does not change the Zn II column density, indicating as observed in other GRB 
afterglows (e.g. Savaglio & Fall 2004) that the contamination of the corresponding Zn II 
line by Mg I is negligible. With a Zn II column density of 10 14 3±a3 cm" 2 , GRB 050401 also 
has the highest column density Zn II absorber ever recorded in a DLA system — compare 
to io 13 - 95±0 05 cm" 2 in the previous record holder, GRB 990123 (Savaglio et al. 2003), or 



1 See http://www.ucolick.org/~xavier/SDSSDLA/tab^llfit.html 
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10 i3.4±o.i cm -2 in q S0 1850+40 (Prochaska & Wolfe 1998), the highest value so far detected 
in a QSO-DLA. For GRB 050401 we find [Zn/H]= -1.0 ± 0.4 (taking into account the 0.3 
dex uncertainty on the determination of logiV(Hi)). This metallicity is comparable to the 
few previous measurements for GRBs at similar redshifts, but much higher than those of 
the strongest iV(H I) GRB-DLAs observed to date (e.g. Savaglio et al. 2003; Vreeswijk et al. 
2004; Prochaska et al. 2004; Starling et al. 2005; Chen et al. 2005). This optically determined 
metallicity is lower than the metallicity obtained from the X-ray column density. The best-fit 
equivalent hydrogen column density found from the X-ray data is logiVn = 22.21+g;og for a 
Solar metallicity model, which implies a metallicity lower than Solar ([X/H] = —0.4 ± 0.3) 
when compared to the neutral hydrogen column density, logiV(Hl) = 22.6 ±0.3, determined 
above. This is discussed below in § 4. 
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Table 2: List of all absorption lines detected with S/N>3 



Observed 


Observed 


Probable 


Wavelength (A) 


Equivalent Width (A) 


ID 


5074.4 


13 ±3 


I 1302; z = 2.8992 


5178.3 


7 ± 2 




5194.7 


7 ± 2 




5206.0 


10 ±3 


C II 1334, C II* 1335; z = 2.8992 


5318.2 


6 ± 2 




5422.2 


7 ± 2 


C IV 1548; z = 2.4972 


5431.9 


7± 2 


C IV 1550; z = 2.4972 Si IV 1393; z = 2.8992 


5841.5 


9 ± 2 


Al II 1670; z = 2.4972 


5900.0 


6 ± 2 




5950.4 


9 ± 1 


Si II 1526; z = 2.8992 


6036.5 


5 ± 1 


C IV 1548; z = 2.8992 


6044.6 


5 ± 1 


C IV 1550; z = 2.8992 


6273.9 


6± 1 


Fe II 1608; z = 2.8992 


6559.9 


4 ± 1 




6507.9 


6 ± 1 




6517.6 


4.2 ± 0.8 


Al II 1670; z = 2.8992 


<—rr\ a r\ r\ 

7049.0 


3.2 ± 0.8 


Si II 1808; z = 2.8992 


7229.5 


6.0 ± 0.6 


Al III 1854; z = 2.8992 


7262.0 


5.6 ±0.8 


Al III 1862; z = 2.8992 


^-rn r\r\ c\ 

7899.2 


5.5 ± 0.9 


Zn II 2026; z = 2.8992 


8022.6 


3.4 ± 0.9 


Cr II 2056; z = 2.8992 


8042.2 


4.9 ± 0.7 


Zn II 2062, Cr II 2062; z = 2.8992 


8061.6 


2.1 ±0.7 


Cr II 2066; z = 2.8992 


8197.9 


5.4 ±0.7 


Fe II 2344; z = 2.4972 


8303.8 


4± 1 


Fe II 2374; z = 2.4972 


8331.5 


6.4 ±0.6 


Fe II 2382; z = 2.4972 


8815.8 


3.7 ±0.9 


Fe II 2260; z = 2.8992 


9023.8 


3.7 ±0.6 




9093.7 


9± 1 


Fe II 2600; z = 2.4972 


9137.6 


6.2 ±0.8 


Fe II 2344; z = 2.8992 
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Table 3: Absorption line fits at z = 2.8992 

Central (A) log TV ID [X/H] [X/Zn] 



Wavelength 


(cm 2 ) 








7049.8(6) 


16.5(4) 


Si II 1808 


-1.7(5) 


-0.7(5) 


7900.3(7) 


14.3(3) 


Zn II 2026 


-1.0(4) 




8042.7(7) 


14.3(3) 


Zn II 2062 






6271.7(5) 


16.0(2) 


Fe II 1608 


-2.1(4) 


-1.1(4) 


6282.4(5) 


16.0(2) 


Fell 1611 






8772.7(8) 


16.0(2) 


Fe II 2249 






8815.2(8) 


16.0(2) 


Fe II 2260 






8017.8(7) 


14.6(2) 


Cr II 2056 


-1.7(4) 


-0.7(4) 


8041.1(7) 


14.6(2) 


Cr II 2062 






8056.4(7) 


14.6(2) 


Cr II 2066 







Note. - 



The uncertainty in the last significant digit is listed in parentheses after that digit. 
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Fig. 5. — Simultaneous Voigt-profile fitting of carefully selected, weaker and unblended lines 
observed at z^ s = 2.8992 in the afterglow spectrum of GRB 050401. We also show two more 
saturated/blended lines (Si II A1526 and Fe II A2344) that were not used in the line fitting. 
The origin of the velocity scale corresponds to the redshift z a b s = 2.8992. 
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We note that there is an intervening absorption system at z — 2.4972 based on Si II, 
Al II, and Fe II lines. The metal lines in this system are strong but significantly weaker than 
those from the GRB system. Moreover, the Zn II lines from this system are not detected. 
It is not possible to estimate the H I column density from the Lyct line as there is no flux 
left from the afterglow around 4255 A. Assuming a metallicity of a tenth solar we estimate 
logiV(Hl) ~ 21. This intervening DLA therefore provides a negligible contribution to the 
uncertainties on the absorbing column in X-rays or the dust extinction. 

In the 2 dimensional spectrum we also identify a Lya emitter at z — 2.65 at an impact 
parameter of 10.2" south of the GRB position (81.3 kpc at z — 2.65, see the upper panel in 
Fig. 1). There is no significant absorption in the afterglow spectrum at this redshift. There 
is also no evidence for Lyct emission from the host galaxy in the 2D spectrum, but as the 
host galaxy is not detected this does not exclude even a large equivalent width Lya emission 
line from the host. 



3.3. The spectral energy distribution 

The photon-weighted mean time of the WT-mode X-ray spectrum was 1075 s after the 
burst. The flux of the optical spectrum was interpolated to this time, using the cto = 0.82 
optical decay slope and the flux calibration of the spectrum. The best-fit to the R-band 
lightcurve is consistent with the flux interpolation of the optical spectrum to this time 
(Fig. 4). The NIR data were also corrected to this time, assuming their decay was the same 
as the R-band decay rate. The correction in the NIR was less than for the optical spectrum 
since the NIR observations occurred between 37 and 76 minutes after the burst. 

The SED (Fig. 4) clearly shows that a break in the optical-to- X-ray spectral index 
is required even to fit a simple SMC extinction-correction to the optical spectrum alone. 
Larger absorption with less reddening would of course allow the optical spectrum to be fit 
if we assume an intrinsic optical slope that was a continuation of the X-ray slope. Such an 
extinction would be consistent with the idea of a flat or 'grey' extinction as advocated by 
Galama & Wijers (2001), Savaglio & Fall (2004) and Stratta et al. (2005) among others. 

Limits on the total extinction (A\) for GRB 050401 were derived from the extreme cases 
of 1) a single power-law extrapolated from the X-ray spectrum to the optical {(3 = 0.89) and 
2) a broken power-law with a A/5 = —0.5 (giving j3 = 0.39), as expected from the fireball 
model (Sari, Piran, & Narayan 1998), with the break set conservatively just inside the 
observed X-ray band (0.4 keV). These limits on A A are shown in the inset to Fig 4. It is 
notable that even a very flat extinction curve cannot alone reconcile the upper limit to the 
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extinction with the large metal column densities we observe. We discuss this conflict below. 

4. Discussion 

The optical afterglow of GRB 050401 is particularly faint — in fact it is the faintest GRB 
afterglow for which an absorption line redshift has been obtained. The observed optical to 
X-ray spectral slope is flatter than expected for the fireball model, therefore GRB 050401 
falls into the 'dark burst' population as defined by Jakobsson et al. (2004a). The observed 
reddening in the optical spectrum allow us to conclude that the 'dark' nature of the burst 
can be explained by extinction in the host galaxy. It is interesting that the very high column 
densities observed in GRB hosts, which are probably local to the GRBs, means that we 
do not need to invoke very dusty environments to explain dark bursts — even with a small 
dust-to-gas ratio these very high column densities are sufficient to darken the afterglows. 

4.1. Metallicity 

Since we expect GRB progenitors to be short-lived (Hjorth et al. 2003b; Stanek et al. 
2003), the cloud probed by the afterglow is likely to be the matrix out of which the progenitor 
was born. For this reason, it may be problematic to try to infer the global properties of 
the ISM at high-redshift from the GRB's probable host DLA. But for exactly the same 
reason GRB-DLAs are particularly important for determining the nature of the progenitor 
material of the GRB. For instance, the favoured model of long duration GRBs (the 'collapsar', 
Woosley 1993; Woosley & Heger 2006), predicts that GRBs should generally be found in low 
metallicity environments. Alternative progenitor models exist where there is no such low 
metallicity preference (e.g. Ouyed et al. 2005; Fryer & Heger 2005). 

Previous observations of GRB afterglows and GRB host galaxies have indeed suggested 
that GRBs reside in low metallicity galaxies (Le Floc'h et al. 2003; Fynbo et al. 2003; 
Prochaska et al. 2004), and in fact the two next strongest GRB-DLAs after GRB 050401 have 
rather low metallicities (Vreeswijk et al. 2004; Starling et al. 2005; Chen et al. 2005). The 
absorbers in GRB afterglows usually have high metal column densities, but relatively little 
reddening — this phenomenon has been attributed to a flat extinction law, often suggested to 
be due to destruction of dust grains (Galama & Wijers 2001; Fruchter et al. 2001; Savaglio 
& Fall 2004; Stratta et al. 2005). 

Here we have one of the first direct measures in a single GRB of the H I column density 
and the total metal (primarily a-element) column density. If we assume that the X-ray 
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and optical absorptions probe the same matter, and there is no evidence to the contrary, 
this allows us to constrain the metallicity from the X-ray measurement, suggesting [X/H] = 
— 0.4±0.3, where most of the uncertainty is in the hydrogen column. Using the best-fit to the 
optical absorption lines we find [Zn/H]= —1.0 ± 0.4. Zn is one of the elements least affected 
by depletion onto dust grains, at least in warm clouds in the Galaxy (Savage & Sembach 
1996), and is also a good tracer of the Fe-group metallicity. By fixing the hydrogen column 
density at the level of the optically-determined H I column and the Fe-group elements at the 
level determined from the Zn abundance, we can fit for the a-element abundance in the X-ray 
spectrum. We obtain a = 0.36 ± 0.06 Z Q) which translates to [a/H]= —0.4 ± 0.3, including 
the H I uncertainty. We can also directly obtain the ratio of a-chain metallicity (from the 
X-ray spectrum) to Fe-group metallicity (from the optical spectrum). This ratio may hint 
at a possible overabundance of a-elements to Zn compared to solar: [a/Zn]= 0.6 ± 0.3 (a 
factor of about 4-5). 

In an analysis of QSO absorption lines systems, Turnshek et al. (2003) invoke an over- 
abundance of a-chain to Fe-group elements of 2.5. Such an overabundance is consistent 
with the results obtained here. Indeed, this degree of a-element overabundance is consistent 
with the Zn metallicity (see Turnshek et al. 2003; Matteucci & Chiappini 2005), and would 
indicate a relatively young ISM still dominated by products from core collapse rather than 
type la supernovae (Matteucci & Chiappini 2005). 

From the optical measurements, there is no evidence for an overabundance of Si, which 
is also considered an a-chain element, compared to Zn, with [Si/Zn]= —0.7 ± 0.5. Since 
Si is depleted onto dust grains as discussed above, this ratio probably reflects differential 
dust depletion effects. Evidence for dust depletion comes from iron for which we measure 
[Fe/Zn]= -1.1 ± 0.4 similar to the highest seen in QSO-DLAs (Ledoux et al. 2003). 

In cold Galactic clouds Zn is depleted by about —0.7 dex (Savage & Sembach 1996), so 
it is important to examine if the apparent alpha-element overabundance could be due to a 
large depletion of Zn onto dust grains. We believe the depletion of Zn is not large in this 
the abundances of Si, Cr, and Fe relative to Zn in this system ([Si/Zn]= — 0.7 ± 0.5, 
[Cr/Zn]= —0.7 ± 0.4, [Fe/Zn]= —1.1 ± 0.4) are more similar to the element ratios in warm 
Galactic clouds ([Si/Zn]~ —0.5, [Cr/Zn]~ —1.0, [Fe/Zn]~ —1-2) than those in cold Galactic 
clouds ([Si/Zn]~ —0.6, [Cr/Zn]~ —1.6, [Fe/Zn]~ —1.6) (see Savage & Sembach 1996, for 
details). Nevertheless, we cannot exclude some depletion of Zn onto dust and in this case 
the total metallicity would be higher if metal ratios were similar to the solar values ([X/H] = 
—0.5 ± 0.3). Strictly, the derived ct-element overabundance should hence be considered an 
upper limit. 



-23- 



4.2. Extinction 

The optical spectrum is well fit using an SMC extinction curve with Ay = 0.62 ± 0.06 
and an intrinsic slope (before dust extinction) of /3 — 0.5 ± 0.2. With a slightly flatter slope 
of j3 — 0.39, as expected from the X-ray slope j3 x = 0.89 and a cooling break of A/3 = 0.5 
between the X-rays and the optical, the required extinction is slightly higher, Ay = 0.67 
(see Fig. 4). 

This reddening-derived estimate is somewhat less than might be expected from a simple 
translation of the optical H I or Zn II columns to extinction. These would respectively 
suggest Ay = 2.3l};2 assuming an SMC dust-to-hydrogen ratio (Bouchet et al. 1985), or 
Ay = 2.2+^, following Savaglio & Fall (2004). The similarity of these latter two estimates 
is not surprising of course since the estimated metallicity in the SMC is similar to the [Zn/H] 
observed here. However, these estimates are strongly dependent on the metallicity and if 
some Zn is depleted onto dust, the metallicity-estimated Ay could be substantially greater. 

A more robust, direct comparison can be made between the optical spectrum and the 
optical-to- X-ray SED, thereby placing bounds on the total and selective extinctions (see the 
inset in Fig. 4). Ay < 0.5 is excluded: it would require a non-power-law continuum in the 
optical which was also substantially below the continuum expected from the X-ray spectrum, 
even with a steep spectral break. On the other hand, we can also exclude Ay > 4.5 since 
it requires not only an extremely flat extinction curve, but also that the optical and X- 
ray afterglows are driven by different emission components (e.g. that the X-ray emission is 
dominated by inverse compton radiation rather than synchrotron radiation). 

The absolute limits on the extinction (4.5 > Ay > 0.5) and the extinction estimates 
based on metallicity and on reddening-fitting examined above, are all much less than the 
Ay inferred from the metal column derived from the soft X-ray absorption, which is 9.1^};g 
(Predehl & Schmitt 1995). Since the X-ray absorption is dominated by metals, the conversion 
from X-ray absorption to dust column is not sensitive to a gas-to-metals conversion, but to 
the dust-to-metals conversion. Allowing an a-element overabundance as suggested above 
would not make these more compatible since it is the a-elements that compose most the 
dust. We are therefore forced to invoke a dust-to-metals ratio (as distinct from the dust-to- 
gas ratio or the metallicity) which is smaller in this DLA than in the Galaxy or the Magellanic 
clouds. If the reddening-derived extinction measure is correct, this is more than a factor of 
ten less than the SMC dust-to-metals ratio. But such a resolution to the low-reddening vs. 
high metal column problem clashes with the apparent pattern of dust depletion observed 
above, as depletion at this level naturally implies the presence of a large dust column along 
with the gas. We are left with a paradox: we very clearly observe a large metal column 
density and insufficient total extinction (and reddening) implying a low dust-to-metals ratio, 
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but at the same time we have some evidence of dust depletion effects in the ratios of the 
optical lines. This is puzzling. We cannot exclude the hypothesis that there is a large dust 
column which has very little effect on the transmission of optical light. On the other hand 
it seems reasonable to say that the evidence for dust depletion (4-5cr in total) may not be 
quite strong enough to support such a far-reaching conclusion. In that case we invoke a low 
dust-to-metals ratio as above. Such a low dust content with a moderately high metallicity 
would indicate either sublimation of most of the dust along the line of sight (e.g. Draine & 
Hao 2002), or more simply, that the dust observed in this DLA system has been formed by 
SNe rather than AGB stars. This is because SNe contribute metals and dust very quickly to 
the ISM, but while they provide nearly all of the metals, they produce only relatively small 
amounts of dust (SNe produce at most only a few percent of the dust found in the local 
universe; Ferrarotti & Gail 2005; Nittler 2005). 



5. Conclusions 

We have reported the strongest damped Lya line observed to date, in the spectrum of 
GRB 050401. The metallicity is fairly high even for a GRB-DLA, with [Zn/H] = -1.0 ±0.4. 
Large X-ray absorption is also detected, from which a large a-element column density is 
inferred. The extinction expected based on the X-ray absorbing column density exceeds 
the total extinction allowed by the SED in GRB 050401, and is more than an order of 
magnitude greater than the extinction inferred from a good fit of the SMC extinction curve 
to the reddening in the optical spectrum. A grey or flat extinction law cannot explain this 
discrepancy. We suggest that since the metallicity is high, and the total extinction is low, it 
implies that there is comparatively little dust along the line of sight to this GRB (however 
there is puzzling evidence to the contrary in the depletion patterns of the optical absorption 
lines). This high metal column with little extinction may be evidence of dust destruction, 
or that the stellar population is very young and that little dust has been formed yet, in 
spite of the vigorous star-formation activity. Based on these results, inferences on dust 
extinction properties drawn from comparisons of optical reddening and metallicity should 
be approached cautiously. 

We thank Peter Capak for excellent support of our CTIO observations. The Dark 
Cosmology Centre is funded by the DNRF. We acknowledge benefits from collaboration 
within the EU FP5 Research Training Network, 'Gamma- Ray Bursts: An Enigma and a 
Tool'. 

We note that De Pasquale et al. (2006) have published an analysis of the afterglow of 
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GRB 050401 while this paper was in review. The analysis of the XRT data by De Pasquale 
et al. produced compatible results with those obtained here. The optical afterglow lightcurve 
published here is more complete and yields a steeper single power-law decay. We also show 
here that grey dust is not a sufficient explanation for the discrepancy between the apparent 
extinction in the optical/UV and the large X-ray column. 
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